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^ ; ABSTRACT 

>, 

OO , The photometric data for a short-period (0.1984 day) eclipsing binary V38 

fvj I discovered by the OGLE micro-lensing team in Baade's Window field BW3 have been 

CN I analyzed. The de-reddened color (V — Ic)o = 2.3 and the light-curve synthesis solution 

1^^ ' of the /-filter light curve suggest a pair of strongly-distorted M-dwarfs, with parameters 

0^ ■ between those of YY Gem and CM Dra, revolving on a tightest known orbit among 

,^ ■ binaries consisting of Main Sequence stars. The primary, more massive and hotter, 

II . component maybe filling its Roche lobe. The very small amount of angular momentum 

o 



in the orbital motion makes the system particularly important for studies of angular 
'^ I momentum loss at the faint end of the Main Sequence. Spectroscopic observations 

^ ' of the orbital radial velocity variations as well as of activity indicators are urgently 

needed for a better understanding of the angular-momentum and internal-structure 
^ ' evolutionary state of the system. 

1. INTRODUCTION 

The authors of this paper share a strong conviction that light-curve solutions of individual 
eclipsing binary stars should belong to the past, with the exception of extraordinarily important 
systems. We feel that the variable star V38 in the third (BW3) field of the OGLE micro-lensing 
project is such a system. With the orbital period of only 0.1984 day, it is the shortest-period 
known binary consisting of Main Sequence stars. 
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BW3 V38, from now on called #3.038 (02000 = 18''04'"44!19, (52ooo = -30°09'05'.'l, 
P = 0.19839 day, Imax = 15.83, {V — Ic)max = 2.45, A/ = 0.78), appears in the second instalment 
of the OGLE Catalog of periodically variable stars (Udalski et al. 1995) among contact binaries 
of the W UMa-type (EW). However, it did not pass the Fourier filter and was not included in 
the impersonally-selected sample of contact binaries in the discussion of all contact systems in 
Baade's Window (the R-sample) by Rucinski (1997 = R97). But this rejection was a marginal one 
because the Fourier filter had been constructed for the most common G-type systems observed in 
the F-band, whereas in this case, it was applied to a M-type system observed in the /-band. It 
can be argued that light curves of contact binary systems are primarily dependent on the strongly 
perturbed geometry of the stars, and weakly dependent on the atmospheric properties (Rucinski 
1993), but it is possible that application of the Fourier filter was carried too far for the case of 
#3.038. A full light-curve solution seemed to be in order. 

The importance of #3.038 stems from its very short period, but also, if it is a contact system, 
from its very late spectral type. Short-period, late-type systems are common among contact 
binaries, but the period distribution of these systems shows a sharp cutoff, which is well defined by 
the system CC Comae, with the orbital period is 0.221 day and the intrinsic color V — Ic = 1.39 
(Bradstreet 1985). In the volume-limited sample to 3 kpc formed from the OGLE sample (R97), 
contact binaries appear with high apparent frequency only within intervals of orbital periods and 
intrinsic colors 0.25 < P < 0.5 day and 0.4 < y — / < 1.2. This color range translates into the 
range of spectral types of about F2 to K5. It is an interesting puzzle of the stellar structure theory 
why the short period cutoff is so well defined and sharp, and why there are no contact binaries 
consisting of late-K and M dwarfs. Some arguments have been presented in favor of the stellar 
structures converging to the fully-convective state (Rucinski 1992), but - quite possibly - the 
reasons for non-existence of very low-mass systems are different and more complex. The discovery 
of #3.038 possibly confronts us with a contact system beyond the current period cutoff. The 
spectral type of #3.038 must be much later than that of CC Com, as the observed colors of the 
two systems are V — Ic = 2.45 and 1.39, respectively. A simple estimate of the reddening based 



on the maps of the reddening of Stanek (1996) (see Sec. |2.lD suggests Ey-i = 0.18. Thus, #3.038 



is much cooler than CC Comae and consists of M-type dwarfs. 

Even if #3.038 is not a contact system, it is extremely interesting object as a very close 
binary consisting of M-type dwarfs. Most stars in the Galaxy are M-type dwarfs, yet ~ because 
of their faintness - we know very little about them. In fact, the low end of the Main Sequence is 
calibrated by only two eclipsing systems, YY Gem and CM Dra. Because of large numbers of the 
M-type dwarfs, small systematic errors made in analyses of these two stars may have profound 
effects in our understanding of even so remotely related subjects as dynamics of the Galaxy or the 
critical density of the Universe. In addition, analysis of old, late-type dwarfs can give us such basic 
data as the primordial helium and metal abundances (Paczyiiski & Sienkiewicz 1984; Chabrier 
& Baraffe 1995; Metcalfe et al. 1996). We should note that CM Dra is a Population II object 
while YY Gem belongs to a sextuple a Gemini (Castor) system which must be relatively young 



(age about 500 Myr) as it contains an unevolved AlV star. Thus, we urgently need more systems 
to clearly understand effects of age and metallicity in the lower Main Sequence. 

The strong distortion of components is directly visible in the light curve (Figure Q) of ^^3.038. 
With continuing angular-momentum loss, which must operate in such a late-type system, the 
binary is apparently on its way to become a contact system. The major question is: Why we do 
not see such systems in large numbers? Is the angular momentum loss so severe that contact 
systems very quickly merge and form single stars? Or, at the other extreme, is the loss so feeble 
that most among close M-dwarf pairs did not lose enough angular momentum in the lifetime of 
the Galaxy to form contact systems? If the latter is the case, is #3.038 one of the oldest among 
such close M-type binaries? 

This paper is meant as a first, exploratory attempt to extract as much information from the 
extant data as possible. We hope to stimulate interest in this important system; in particular we 
hope that the much needed spectroscopic observations will be obtained soon. 



2. THE LIGHT CURVE SOLUTION 

2.1. Assumptions 

While attempting a more detailed study of #3.038, one is confronted with the limited extent 
of the available data. The OGLE data published in the catalogs and available over the ftp service 
consist of the light curve in the I-band, the maximum-light V — Ic color and the orbital period. 
The light curve (Figure ||) is rather noisy, with a standard error at maxima of ~ 0'".02, but with 
deep eclipses (A/ = 0.78) suggesting a high value of the inclination. An analysis of the light 
curve is still sensible, but one has to be careful about the strategy to apply in the solution and 
about the estimation of uncertainties of the derived parameters. With any standard light-curve 
synthesis method, such as the Wilson and Devinney (1971 = WD) program (we used its updated 
1992 version), the main and the most obvious concern would be the choice and the number of free 
parameters. We address this matter below, but first we look into the available information about 
components of #3.038. 

The observed, possibly strongly-reddened V — Ic color is the only information about the 
spectral type of the components. To derive the unreddened V — Ic color and the approximate 
absolute magnitude of the system, we used a simple iterative procedure which followed the one 
presented in R97. The R97 computations were based on Stanek (1996) reddening maps of Baade's 
window and on the calibration period-color-magnitude established for W UMa binaries. We do not 
know whether #3.038 is actually a contact system. Therefore, we used the same iterative approach 
as in R97, but replaced the W UMa calibration with the Main Sequence My = MyiV — Ic) 
relation of Reid and Majewski (1993), shifted upward by 0.75 mag for a binary system of identical 
components. The final values are: {V — Ic)o = 2.3 it 0.1, assuming errors of 0.03 and 0.5 in V — Ic 



and M/, respectively. The predicted systemic absolute magnitude is M/ ~ 7.5 and the distance is 
d = 400 lb 85 pc. The relatively red intrinsic color corresponds to a M3 spectral type, so that a 
reasonable value of primary effective temperature is Ti = 3500 K (Bessell 1979 , 1990). 

Extensive tables of limb darkening coefficients for R, I, J, K bands have been recently 
published by Claret et al. (1995). We used the (linear) limb darkening coefficient, x = 0.56, 
corresponding to the primary effective temperature and log 5 between 4 and 5. The gravity 
darkening and reflection coefficient were fixed at the standard values for convective envelopes 
{g = 0.32, A = 0.5; Lucy 1967, Rucinski 1969). The low effective temperature of primary does not 
allow the use of the model-atmosphere option in the WD code which is limited to the temperature 
range 4,000 - 25,000 K (we note that this option in the 1992 version of the code uses the obsolete 
model atmospheres from 1969). All computations were therefore performed in the black body 
approximation. For a one band-pass light-curve in /, this is definitely not a critical assumption. 

The maxima of the light curve show slightly different heights. The difference is about 
51 = 0.02 which can be modeled by a dark spot on the larger star. In order not to introduce too 
many free parameters into the geometric-elements solution, the spot parameters were determined 
by trial an error before starting the differential correction procedure. Since the spot characteristics 
are very weakly constrained by light curve solutions without spectroscopic information (Maceroni 
and Van't Veer 1993), we assumed the spot to be located on the equator of the primary component. 
The other fixed parameters were: longitude and angular radius of the spot 290 and 11 degrees, 
respectively, and the temperature factor Tspot/Tstur = 0.8. The advantage of having a spot is in 
preventing occurrence of oscillations in the iterative procedure between solutions which fit one 
maximum at a time. 



2.2. Results of the solution 

The input model for the light curve solution consisted of a system of two identical stars with 
the characteristics described in the previous section and exactly filling the inner critical Roche 
surfaces. From this point, the surface potentials were left the freedom to vary in an independent 
way, adopting the WD "mode 2" of operation, which is normal for detached systems. We used 
seven adjustable parameters: the orbital inclination i (in degrees), the mass ratio of the system q 
(expressed as the mass of the eclipsing component to the eclipsed one in the deeper minimum), 
the two surface equipotentials i7i and Q2 (defined as in the WD code), the secondary effective 
temperature T2, the primary component luminosity Li in the spectral band used and the phase 
shift of the primary eclipse Acj) (expressed in units of the orbital period). Their subdivision in 
two subsets (Wilson and Biermann 1976) helped to control the strong inter-dependencies present 
among some of the parameters during the iteration process. The procedure quickly converged to 
the solution presented in Table ||. This table contains the main solution as well as the results of an 
"bootstrap re-sampling" experiment described below. The geometry of the final model is presented 
in Figure |2|. The two components are very close to, but still inside, their respective Roche critical 



surfaces. However, for the more massive and hotter component, the results are consistent at the 
one-sigma level, with the star actually filling its critical inner potential, popularly known as the 
"Roche lobe". 



2.3. Uncertainties of parameters 

It is well known in the community of users of the Wilson-Devinney program that this fine 
code provides excellent light curves, but that the estimates of errors of the adjustable parameters 
are unrealistically small. The reason is partly the strong correlation between the relatively many 
parameters, partly the non-normal distribution of measurement errors. The WD code provides 
the "probable" errors, e,, which are derived by the differential correction routine and are related 
to the standard errors of the linearized least-squares algorithm through ej = 0.6745(Jj. The errors 
(Tj can be used as a measure of the uncertainties only for normal distributions of the errors. 

As explained in the book "Numerical Recipes", Sec.15.6 (Press et al. 1992), the technique of 
bootstrap re-sampling is probably the most useful for estimation of confidence levels for complex 
least-squares solutions. It uses input data re-sampled with repetitions, establishing confidence 
levels from the spread in resulting parameters. However, its application was not possible with the 
current WD code, which ~ by intent - was designed on the principle of "interactive branching" 
during the differential correction procedure (Wilson 1985). Combining the interactive branching 
with the bootstrap re-sampling would mean performing thousands of separate solutions "by 
hand". Instead, we used a simplified bootstrap approach in which 10,000 least-squares solutions 
were made within the minimum already established by a single, iterated solution. Thus, only 
one set of light-curve deviations and parameter derivatives were used. This, by necessity, would 
give us under-estimates of the real errors. However, as we have found through application of 
this approach, these estimates are larger that the formal mean standard errors and do show 
strong inter-parametric correlations. We should remember that inter-parametric correlations are 
not always a detrimental circumstance, as some derived parameters (such as radii) may show 
smaller uncertainties than for the case of simple quadrature addition of uncorrelated errors. The 
bootstrap-estimated errors tend to be also somewhat pessimistic as some re-samplings give light 
curves so clumpy and unevenly covered that nobody would ever attempt solving them. Since we 
see potential for under- and over-estimation of the errors, we give in Table |l| both sets of error 
estimates, the mean standard errors for the iterated solution and the one-sigma ranges around the 
median values for the bootstrap experiment. 

A selection of contour plots giving two-dimension representations of the one-sigma uncertainty 
levels is shown in Figures |3| and §. They have been obtained by finding contours enclosing 68.3% 
of the bootstrap solutions around their respective median values. They are frequently elliptical 
(indicating inter-parametric correlations), but symmetric around the medians of the respective 
seven primary parameters. The symmetry is due to use of one set of derivative within the 
minimum of the summed squared deviations. It should be noted that the bootstrap solutions 



tend to cluster slightly away from the best, iterated solution. This is caused by the fact that the 
bootstrap solutions are always based on light curves with poorer coverage than the full solution. 
This biases the distribution of the mass-ratio q, as this essential parameter is determined solely 
from the distortion of the light curve between minima. If this part of the light curve is poorly 
covered, or if observations have large errors, then the mass-ratio would tend to take smaller values. 
The only proper resolution of this difficulty would be to obtain spectroscopic radial- velocity data 
which would stabilize the solution. 



2.4. Relation to the Roche geometry 

The Figure § contains uncertainty plots for the relative radii of the components. Here we 
used the "side" radii in the orbital plane which - as can be shown - are practically equal to the 
volume radii. In addition to the mean radii from the iterated solution, we also show the median 
values. Note that they are not exactly in the centers of the one-sigma contours, as for the other, 
primary elements. This is due to the inter-parametric correlations, since the radii are given by 
expressions of the type r, = ri{Qi, q). 

As we can see in Figure Q, the results for the radii are not well constrained. However, a plot 
in the same figure of the uncertainty contours for the surface potential of the larger component, 
ill, indicates that the primary almost fills its Roche lobe, as the potential is slightly larger than 
the critical one. The one-sigma contour follows the critical potential curve for the relatively large 
range of mass-ratios permitted by the current solution. If the primary actually fills its lobe, any 
possible mass-transfer between components would be undetectable in the current, single-epoch, 
photometric data. 



3. DISCUSSION OF #3.038 PROPERTIES 

3.1. Comparison to YY Gem and CM Dra 

With the photometric solution of the light curve given in Table |l] we still have rather modest 
information about the system. Four parameters result from the solution (i, q and f^i, ^^2 or ri, 
r2), two describe the relative luminosity characteristics of the components (T2 and Li) and two 
more come from the original observations (V — Ic, P)- On the basis of these, we can attempt to 
place #3.038 relative to other eclipsing M-dwarfs on the Main Sequence. There are only two such 
systems: YY Gem (Castor C), analyzed by Bopp (1974) and by Leung & Schneider (1978) and 
CM Dra analyzed by Lacy (1977) and by Metcalfe et al.(1996). Both are well detached systems 
with orbital periods 0.81 day and 1.27 day, and both show only moderate proximity effects, in 
that no mutual illumination of the components is visible and the stars are practically undistorted 
by tidal forces. In this respect, #3.038 is very different as this is the first known M-type binary 



showing strong distortion of components, directly linked to the very short orbital period of 0.198 
day. The ages of the reference stars are very different: YY Gem belongs to the Castor (a Gem) 
system whose age must be young as it contains hot, unevolved star^, whereas CM Dra is a 
Population II object. Nothing is known at this moment about population characteristics of 
#3.038. 

The color of #3.038, {V — Ic)o = 2.3±0.1, gives us information about its effective temperature 
relative to YY Gem and CM Dra. Both reference systems were discussed by Monet et al. (1992) in 
relation to other intrinsically faint stars. In particular, their Kron-Cousins colors were plotted in 
Figure 11 of Monet et al. which relates the colors to the effective temperature. These unpublished 
data had come from observations at the U.S. Naval Observatory^: YY Gem, V — Ic = 1.92; 
CM Dra, V — Iq = 2.94. Bessell (1990) previously measured CM Dra in the same system and 
obtained V — Ic = 2.92. Thus, judging by the colors, #3.038 is placed about half way in the 
effective temperature between YY Gem and CM Dra. 

An independent evaluation of the relation of #3.038 to YY Gem and CM Dra can be obtained 
from mean densities of the stars. The published masses and radii for YY Gem and CM Dra permit 
direct determination of the component densities in both systems. They are: 3.0 ± 0.4 and 4.1 ± 0.6 
g/cm^ for YY Gem and 21 it 3 and 22 it 3 g/cm^ for CM Dra. For #3.038 we can use the Kepler's 
law, as re-written by Mochnacki (1981): pi = 0.079/(Vi(l + q)P'^) and p2 = 0.079^/(^2(1 + q)P'^) 
g/cm^, with the period P in days. This form emphasizes the orbital period which is known 
practically with no error. The relative volumes of components, Vi and V2 (in orbital units), are 
obviously poorly known as they contain uncertainties in radii in third power. Inserting our results 
for #3.038, we obtain pi = 4.8 it 0.4 and p2 = 6.0 it 3.4 g/cm^. These are evaluations based on the 
assumption that the errors in q and Ri are uncorrelated. The bootstrap sampling estimates which 
explicitly take into account the correlations are less pessimistic: pi = S.IIq'^ and p2 = 5.9l]^'5 
g/cm^. 

Instead of using the derived geometric parameters, we can note that the components of 
#3.038 are apparently just inside their critical Roche lobes. Therefore, relatively firm lower limits 
on the densities can be estimated by using the volumes of these lobes, Vi{q) and V2{q). The 
only source of uncertainty in these estimates is then in the mass-ratio q. We used the formulae 
of Eggeleton (1983) to determine the Roche-lobe volumes, Vi{q). Using the iterated solution, we 
obtained pi > 4.17ib0.07 and p2 ^ 4.60±0.05 g/cm^, whereas the booststrap sampling experiment 
gave: pi > ^m^oir and p2 > 4.6611^:1? g/cm^. 

The results on the density clearly show that the densities of the components of #3.038 are 



^According to probably the best extant spectral classification of Garrison & Beattie (1997), the hot components of 
the systems have spectral types Aim, A2Va, A2m, A5V:. A simple comparison with the Pleiades (where the earliest 
spectral types are B6IV and B8V) and Hyades (where the earliest spectral types are A2m and A5) with the respective 
ages of about 80 and 600 Myr, indicates that the age of sextuple system of Castor must be roughly 500 Myr. 

''We are indebted to Dr. Conard Dahn for sending us the data. 



larger than those of the YY Gem and that the primary component in the system may have 
relatively lower density than its companion, which which would be consistent with it being slightly 
more evolved of the two. However, it is doubtful if the densities of the stars in #3.038 are as 
high as those in CM Dra because the firm lower limits are probably very close to the actual 
values. Thus, we have another indication that in terms of the location on the Main Sequence, the 
components of 7^3.038 are in between those for YY Gem and CM Dra. 



3.2. Angular momentum of #3.038 

The short period of #3.038 brings up the question of the angular momentum (AM) of the 
system, H. Disregarding the spin angular momenta which an contribute about 10% to the total AM, 
the formula for the orbital angular momentum is if = 1.2 x 10^^g/(l + g)^(P/l dy/^{Mtot/MQ)^/^ 
in cgs units. For the current period of 0.2 day, and assuming Mtot = 1 Af©, the systems contains 
now about H ~ 1.8 x 10^^ cgs. This a very small amount, obviously due to the small masses and 
the short period. The question is: How much AM has been lost by the system in its evolution? 
The answer depends on when, at which orbital period, the components became inter-locked in 
the spin-orbit synchronism. If we assume that this happened at about P ~ 2 days, which seems 
a plausible number, i.e. when H ~ 3.9 x 10^^ cgs, and simply subtract the two values, then the 
answer is that #3.038 has lost as much AM in evolving from 2 days to 0.2 day as it contains 
now. In order to evaluate the AM loss rate, we should know how long has this process has taken. 
For this we have no answer as we have no idea about the age and evolutionary state of the 
system. However, since the AM loss must have been appreciable, studies specifically addressing 
the questions of age and evolution will definitely shed light on the totally unexplored matter of 
the AM loss for rapidly-rotating M-dwarfs. 

The angular momentum loss (AML) is normally thought to be related to the overall level of 
magnetic activity of stars, but the relation is still poorly known. It appears that both, the activity 
and the AML tend to be "saturated" at high rotation rates (Vilhu 1987). St§pieh (1995) presented 
a formula for the AML in solar-type stars which explicitely relates it to the X-ray activity. 

Results of further studies of #3.038 may have important ramifications ranging from the AM 
evolution soon after reaching the Main Sequence (Hartmann & Noyes 1987, Collier Cameron &: 
Jianke 1994, Stauffer 1994, Soderblom 1996) to the famous period-gap in cataclysmic variable 
stars (Rappaport et al. 1983, Spruit & Ritter 1983, King 1994). In the case of earlier, solar type, 
components the secular evolution of the orbital period, as produced by AML and AM transfer by 
spin-orbit synchronization has been studied by several authors (Vilhu 1982, Maceroni and Van't 
Veer 1991, Maceroni 1992). 

At present, we have no measurements of activity of #3.038. A deep, pointed observation with 



the ROSAT satellite (Guinan, private communication^), shows no X-ray emission at the place of 
the star with the upper limit of about 0.003 counts/sec, which translates into the upper limit to the 
X-ray flux: fx < 1-8 x 10~^^ erg/cm^/s. This upper limit is unfortunately about ten times higher 
than the expected level of X-ray emission from 7^3.038 assuming the 0.1% bolometric-to-X-rays 
conversion efficiency which maybe typical for most active late-type dwarfs (Fig. 17 in Stauffer et 
al. 1994): For Imax = 15.83, assuming the bolometric correction BCi ~ 0.5 for y — / ~ 2.3 (Bessel 
1991), we obtain mtoi — 15.3 and the bolometric observed flux fi,oi — 1.9 x 10"^^ erg/cm^/s. 
Therefore, the expected X-ray flux, using the Stauffer et al. scaling, would be about 1.9 x 10~^^ 
erg/cm^ /s. Additionaly, this estimate does not include the neutral- hydrogen absorption which 
maybe large in this direction, at the distance of about 400 pc. It is known (Pye et al. 1994), 
however, that binary systems - even wide ones ~ can have considerably elevated X-ray activity, so 
that the matter is definitely not closed. 



4. CONCLUSIONS 

We have presented analysis of the /-band light curve and of the photometric colors for the 
closest known pair of M-type dwarfs orbiting each other with the period of less than 0.2 day. 
The components of #3.038 appear to have properties intermediate between those of YY Gem 
and CM Dra; what sets them apart from these well-studied stars, is their strong distortion. The 
photometric solution indicates that the more massive component is very close to or at it critical 
Roche lobe. The present solution is however of moderate quality as the crucial parameter of the 
mass-ratio is poorly constrained by the photometric data. A spectroscopic radial-velocity study of 
this important system is very much needed. 

Tightness of the orbit and the very late spectral types of the components invoke a series of 
questions related to the angular momentum loss in very close M-type dwarfs and to non-existence 
of contact binaries with spectral types later than middle-K. Why we do not see more systems like 
#3.038? Is the angular momentum loss too strong or too weak to produce observable quantities of 
such very close systems? What is the fate of the system with the continuing angular momentum 
loss? Answers to these questions require further observational data for this interesting system. 

We would like to thank Drs. Conard Dahn, Ed Guinan, Kazik St§pieh and Brian Beatty for 
useful comments, suggestions and information. The OGLE team is thanked for the open access to 
their database. 

The research grant from the Natural Sciences and Engineering Council of Canada to SR is 
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^We are indebted to Dr. Ed Guinan for analysis of the Rosat data. 
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Fig. 1. — Light curve of ^^3.038 in the /-band from the OGLE database. The continuous Une gives 
the best fit. The deviations are shown in the panel above the hght curve. 
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Fig. 2. — The three-dimensional picture of ^^3.038 at the orbital phase 0.25. 
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Fig. 3. — Results of 10,000 "bootstrap re-sampling" solutions, as described in the text. The 
contours define the one-sigma uncertainty levels by enclosing 68.3% of all cases. The poor definition 
of some of the contours is entirely due to the contouring routine. The large crosses give the results 
of the mean iterated solution, with lengths of arms equal to the standard errors of parameters. 
The small slanted crosses give the median values for individual parameters. Refer to the text 
for the reasons why the mean and median values do not coincide. Note that the temperature of 
the secondary component was obtained by assumption of the fixed temperature of the primary, 
Ti = 3500 K. 
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Fig. 4. — Additional combinations of parameters in the "bootstrap re-sampling" solutions. The 
broken lines refer to critical (Roche) surface values. Since the radii are the derived parameters, 
their median values are sometimes shifted from the centers of the one-sigma contours. 
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Table 1. Light curve solution of BW3_V38 

Parameter Solution Bootstrap 

mean & std err median & 1-sigma 

i (degr) 85.7 ±1.2 85.8 ± 2.2 

q 0.77 ±0.09 0.70 ±0.18 

^i 3.48 ±0.13 3.40 ±0.32 

^2 3.55 ± 0.25 3.36 ± 0.48 

Ts (Ti = 3500) 3459 ± 15 3456 ± 28 

Li/{Li+L2) 0.597 ±0.036 0.597 ± 0.065 

A0 (phase units) ±0.0020 ± 0.0006 ±0.0018 ± 0.0012 



ri (side) 




0.380 ±0.011 


0.382+-- 


r2 (side) 




0.322 ± 0.050 


n Q9O+0.029 

U.OZO_QQ3g 


pi (g/cm^) 




4.8 ±0.4 


5.1^7 


P2 (g/cm^) 




6.1 ±3.4 


5.9+11 


pi (g/cm^) 


lim. 


> 4.17 ±0.07 


> 4.09+[J-l^ 


P2 (g/cm3) 


lim. 


> 4.60 ± 0.05 


> 4.66l°:l? 



Note. — ri, r2, pi, p2 have been derived from the solution 
of the first 7 parameters. See the text for explanations of the 
lower limits to pi and p2- 



